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ABSTRACT 

The Vela and Geminga pulsars are rotation powered neutron stars, which have been 
identified in various spectral domains, from the near-infrared to hard 7-rays. In the 
near-infrared they exhibit tentative emission excesses, as compared to the optical 
range. To check whether these features are real, we analysed archival mid-infrared 
broadband images obtained with the Spitzer Space Telescope in the 3.6-160 /im range 
and compared them with the data in other spectral domains. In the 3.6 and 5.8 /im 
bands we detected at ~ (4-5)<r significance level a point-like object, that is likely to 
be the counterpart of the Vela pulsar. Its position coincides with the pulsar at < 0'/4 
lcr-accuracy level. Combining the measured fluxes with the available multiwavelength 
spectrum of the pulsar shows a steep flux increase towards the infrared, confirming 
the reality of the near-infrared excess reported early, and, hence, the reality of the 
suggested mid-infrared pulsar identification. Geminga is also identified, but only at 
a marginal 2a detection level in one 3.6 ^m band. This needs a farther confirmation 
by deeper observations, while the estimated flux is also compatible with the near- 
infrared Geminga excess. The detection of the infrared excess is in contrast to the 
Crab pulsar, where it is absent, but is similar to the two magnetars, 4U 0142+61 and 
IE 2259+586, showing similar features. We discuss X-ray irradiated fall-back discs 
around the pulsars, unresolved pulsar nebula structures, and pulsar magnetospheres 
as possible origins of the excesses. We note also possible infrared signatures of an 
extended tail behind Geminga and of the Vela plerion radio lobes. 
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1 INTRODUCTION 

The two nearby middle-aged pulsars, Vela (PSR 
B0833-45) and Geminga (PSR J0633+1746), be- 
long to a small group of rotation powered neutron 
stars (NSs), which are identified in various spectral 
domains outside the radio range. Both pulsars have 
counterparts in the near- infrared (IR) and optical (e.g., 
Shibanov et"aH 120031 , l200fj ; iMignani. Zharikov fc Caraveol 



20071. and references therein ), in th e ultraviolet (UV) 
jRomani. Kargaltsev &. Pavlov 20051: Kargaltsev et al 



20051) , soft X-rays (e.g., iPavlov et al.ll200 j : lKargaltsev eta! 



2005; C araveo et al J 12004, a nd references th erein), hard 



X-rays l|Harding et al l l2002i ) 



and 7-rays (|Abdo et al.l 
l2010al lbl. and references therein). Emission mechanisms of 
pulsars are still poorly understood, and the presence of 
such multiwavelength objects is very useful for the progress 
in this field. 

In soft X-rays and extremal ultraviolet (EUV) the emis- 
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sion from the Vela and Geminga pulsars is dominated by the 
thermal blackbody (BB) like radiation believed to be origi- 
nated from cooling surfaces of the NSs, and, possible, from 
their hot polar caps. In the rest ranges the spectra are de- 
scribed by a power-law (PL) with different spectral indices 
in different spectral domains. This emission is nonthermal 
and generated by relativistic particles accelerated in mag- 
netospheres of the rapidly rotating and strongly magnetized 
NSs. 

In the optical the nonthermal spectra of both pulsars 
are significantly flatter than in soft X-rays, and the optical 
fluxes are considerably below the long wavelength extrap- 
olations of their nonthermal X-ray spectra. This suggests 
a spectral break between the two ranges. Similar spectral 
picture is observed for the young Crab pulsar. In the near- 
IR the Vela and Geminga pulsars exhibit a pparent flux ex- 
cesse s , as c ompared to the optical range ( Shibanov et alj 
I2003L hood ). This is in contrast to the Crab, whose near 
and mid-IR fluxes are compati ble with the spectral extrap- 
olation from the optical range ijSandberg fc Sollermanl l2009; 
iTemim et al.|[200g ) . The near-IR excesses appear to increase 
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Table 1. Spitzer observations of the Vela and Geminga pulsars. 



AOR key 


Dev. 


Obs. date, 


Channels, 


Int. time", 






yy/mm/dd 


fim 


s 


the Vela pulsar 


11374848 6 


IRAC 


04/12/17 


3.6, 5.8, 


72 








4.5, 8.0 


72 


11375104 6 


MIPS 


05/03/07 


24 


92 








70 


335 








160 


105 


11542784 c 


IRAC 


04/12/18 


3.6, 5.8 


9400 






the Geminga 


pulsar 




12540928 d 


MIPS 


05/04/04 


24 


501 


12543232 d 


IRAC 


05/03/26 


3.6, 5.8, 


270 








4.5, 8.0 


270 


19037696 e 


IRAC 


06/10/30 


3.6, 5.8 


14700 


19037952 e 


IRAC 


06/11/22 


3.6, 5.8 


14700 



a The integration time shown for the IRAC is the same for 

the each channel in a given AOR, 

6 PI, Patrick O. Slane, Program ID = 3647 

c PI, Divas Sanwal, Program ID = 3696 

d PI, Michael Jura, Program ID = 23 

e PI, Deepto Chakrabarty, Program ID = 30822 



with the wavelength, and mid-IR observations can be useful 
to confirm or discard their reality. 

Besides the Crab pulsar, there are only two iso- 
lated NSs detected in the mid-IR. Both are anomalous 
X-ray pulsars (AXPs), 4U 0142+61 and IE 2259+586. A 
strong mid-IR flux excess was obse rved for the first one 
ijWang. Chakrabarty fe Kaplan! l200rjh . and si milar, but less 
evide nt feature was found for the second (|Kaplan et al.l 
2009). The excesses can be interpreted as signatures of the 
so-called fall-back discs which can be formed around NSs 
shortly after the supernova (SN) explosions. Such discs can 
be progenitors of pulsar planetary syste ms, one of which 
was discovered around PS R B1257+12 (|Wolszczanl 1 19941 ; 
iKonacki fc Wolszczanll2003l l . However, the association of the 
excesses with the intrinsic pulsar emission still can not be 
excluded. 

In addition to their multiwavelength emis- 
sion, the Vela and Geminga pulsars power ex- 
tended pulsar wind nebu lae (PWNe) most clearly 
in 



200ll; Pavlov et al. 


12001a; Caraveo et all 


2002; 


de Luca et all 12006 


; IPavlov, Sanwal & Zavlinl 


2006: 


Pavlov, Bhattacharyya & Zavlin 2010). The PWN 


mor- 



phologies are different. The younger Vela (~ 10 4 yr old) is 
associated with the supernova remnant (SNR) G263.9— 3.3 
and, as the Crab pulsar, forms a torus-like inner PWN 
with polar jets. The PWN of the older Geminga (~ 
3.4x10 yr) has a tail and/or a bow-shock morphology. 
At larger spatial scales the Vela PWN has a cocoon, or a 
plerion, structure. It is mainly elongated southwards the 
pulsar, that is roughly ortho gonal to the orientation of the 
inner PWN symmet ry axis ( Mar kwardt fc Ogelmanl 1 19951 ; 
iDodson et ail l2003al l. The cocoon has been detected in 
the GeV range with Fermi l|Abdo et all I2010d 1 and the 
AGILE JPellizzoni et al.ll2010l), an d in the TeV range with 
the HESS ( Aharonian et alj " 20061. It is speculated, that 



the high energy source may be a 'relict' PWN distorted by 
the passag e of the remnant reverse shock i n inh omogeneous 
medium (|Blondin. Chevalier fc Friersonl 1200 il l. Geminga 
is not associa ted with any SNR, although, it is discussed 
(|Salvati|[2O10l l, that a dissolved and not visible Geminga 
SNR could be a plausible source of a small scale excess in 
nmlti-TeV cosmic ray p rotons arrival dir ection (spot A) 
discovered with Milagro l|Abdo et alj|2008h . 

Both PWNe have no obvious optical counterparts. 
Some faint extended features were detected in the near- 
IR with the VLT in the l'/ 5-3" 1 vicinity of the Vela pul- 
sar (|Shibanov et al.l I2003T ). They are projected on and/or 
aligned with the south-east jet and the inner arc of its X- 
ray PWN. There is also a faint optical nebulosity with a 
size of few arcsec around the Geminga pulsar detected in 
a deepest up to date optic al image obtained in the I band 
with the Subaru telescope (IShibanov et al.l l2006l. The neb- 
ula is extended perpendicular of the pulsar proper motion, 
and, probably, is the brightest inner part of the bow-shock 
structure which is seen in X-rays. However, the PWN origin 
of those features still needs to be confirmed. 

Motivated by further studies of the IR flux excesses 
in the emission of both pulsars and by searching for their 
PWN IR-counterparts, in this paper we analyse unpublished 
archival mid-IR images of the pulsar fields obtained with the 
Spitzer Space Telescope. We report on finding likely mid-IR 
counterparts of the pulsars and possible features of their 
PWNe. The mid-IR observations are described in Sect. 2, 
the results are presented and discussed in Sect. 3 and 4. 



2 OBSERVATIONS AND DATA ANALYSIS 
2.1 Spitzer data 

The Vela and Geminga pulsar fields have been imaged sev- 
eral times with the Spitzer telescope using the Infrared Ar- 
ray Camera (IRAC) and the Multiband Imaging Photometer 
for Spitzer (MIPS jj. Various instrument channels have been 
used with the effective wavelengths of 3.6, 4.5, 5.8, 8, 24, 70, 
and 160 /im. The information on the respective Astronom- 
ical Observation Requests (AORs) is presented in Table [1] 
We retrieved the data from the Spitzer archive. We used 
the IRAC post-BCD (Best Calibrated Data) mosaic images 
with an effective image scale of 0'.'6 per pixel and the field 
of view (FOV) of 7!5x7?5. They have been generated with 
the Spitzer pipeline procession version of S18.7.0. We also 
regenerated the MIPS post-BCD calibrated mosaic images 
from respective BCD data using the MOPEXBtool vl8.3.1. 
The resulting image scales are 2'.' 4 per pixel for the 24 /im 
MIPS channel with the FOV of 13'3 x 12^8. After the pipeline 
drizzling of multiple individual exposures the effective pixel 
scale for the IRAC mosaics is twice smaller than the nomi- 
nal one (l'.'2 per pixel), while for the MIPS it is consistent 
with the nominal pixel scale of 2 / .'4. 

First inspections of the archival data showed that a 
point-like counterpart candidate of the Vela pulsar is de- 
tected at > 4a significance level in the 3.6 and 5.8 /im im- 
ages of the AOR 11542784, where the integration time was 



1 see |http: / /ssc.spitzer.caltech.edu| 

2 see http://ssc.spitzer.caltech.edu/postbcd/mopex.html 
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maximal among the other AORs for this target (cf. Table 
[l]). The respective mosaic coverage maps show that the real 
integration times of the pulsar neighbourhood located in the 
centre of the mosaic images are only 1-3 % smaller than the 
AOR nominal archival integration time shown in Table[T] At 
the shorter AOR 11374848 the source is probably marginally 
detected at 3.6, 4.5, and 5.8 ^m, and it is not detected at 8 
fim. A bright pixel at the source position is seen in the 24 
/jm image (AOR 11375104), while at longer wavelengths no 
signatures of the pulsar are detected. 

A likely counterpart of the Geminga pulsar is 
marginally detected at ~ 2a significance level in the 3.6 
fj,m images of the long integration time AORs 19037696 and 
19037952. These AORs differ by a zodiacal background level, 
that was ~ 30 % smaller in the second AOR. At 5.8 fim the 
candidate is not resolved even in the image combined from 
the data of the two AORs. However, we found some sign 
of it in the 4.5 /im image of the shorter AOR 12543232. In 
other Spitzer bands and shorter AORs we have not found 
any signature of the pulsar. 

To check that the detected counterpart candidates are 
not the detector artefacts or results of a poor cosmic ray re- 
jection, we analysed the respective badpixel masks, provided 
by the archive together with the data, and regenerated the 
mosaic images from the BCD data using the MOPEX tool 
vl8.3.1. We used the nominal IRAC pixel scale and a half 
of it, as has been done in the pipeline processed images dis- 
cussed above. For the Vela AORs we obtained the similar 
results, confirming the reality of the counterpart candidate 
detection in both IRAC bands of the AOR 11542784 and the 
presence of a bright pixel at the source position in the 24 

image of the AOR 11375104. There is a bright star lo- 
cated in 0.'5 north-east of the Vela pulsar. In 3.6 and 5.8 /im 
images it is oversaturated and produces typical stray light 
artefacts. However, the nearest one is in ~ 8 arcsec south- 
east of the pulsar position and cannot induce a false source 
detection at least in a few arcsec region around the pulsar. 
The marginal Geminga counterpart in 3.6 /im images of the 
AORs 19037696 and 19037952 was found to be present at the 
same significance level in the nominal and the half pixel scale 
regenerated mosaics. However, at the 4.5 /im it disappears 
in the nominal pixel scale mosaic. Its apparent presence in 
the half pixel scale drizzled mosaic is likely to be an artifi- 
cial background fluctuation resulting from an inappropriate 
use of the 0'.'6 pixel scale in the pipeline data procession for 
a shallow observation in this band, where only nine sepa- 
rate exposures have been obtained. Therefore, only an upper 
limit on the Geminga flux can be set in this band. Further 
investigation of the 3.6 /im images of Geminga showed that 
they contain some artefacts. These are so called 'muxbleed' 
and 'muxstripes', which are due to the presence of two bright 
stars in the FOV north-west of Geminga and which cannot 
be perfectly corrected by the reduction tool at given AOR 
setupfl The respective stripes cross the whole mosaics from 
north to south approximately. However, the nearest one is 
located in ~ 0!5 west of the Geminga position and cannot 
result in a false source detection, at least in several tenths 
arcsec pulsar neighbourhoods. 



3 For details, see, e.g., the IRAC Instrument Handbook, 
http://ssc.spitzer.caltech.edu/irac/iracinstrumenthandbook/ 



Thus, we can conclude, that the Vela counterpart can- 
didate in the 3.6 and 5.8 /im images is definitely the real 
source. The Geminga candidate in 3.6 /im images is likely 
to be the real source as well, but it is much fainter and de- 
tected at a lower significance level. The source reality is ad- 
ditionally supported by its independent detection at about 
the same significance level in the two AORs, which were ob- 
tained at different epochs and conditions. Below we mainly 
focus on the data where our counterpart candidates have 
been checked to be the real mid-IR sources. 



2.2 Astrometry 

The coordinates and proper motions of the Vela and 
Geminga pulsars are known with a high accuracy based 



on the radio and optical observations (ICaraveo et alj 
iDodson et al. I [20031; iBignami. Caraveo fc Mereghettil 



2001; 



1993; 



iFahertv. Walter &: Andersonl 120071 ) . For correct positional 
identifications of the pulsars with the sources detected in 
the mid-IR, one has to be sure in the Spitzer pointing accu- 
racy. To check this in the IRAC images of both pulsars, we 
selected about two dozens of unsaturated isolated astromet- 
ric standards from the USNO-B1 catalogue. Their catalogue 
and image positional uncertainties for both coordinates were 

< 0'.'2 and < 0'.'5, respectively. IRAF ccmap/cctran tasks 
were used to find plate solutions. Formal rms uncertainties of 
the astrometric fit were < 0'.'16 with maximal residuals of < 
Cf.'4 for both coordinates. After astrometric transformations 
the shifts between the original and transformed images were 

< ((.'4 for the Vela and < 0'.'2 for the Geminga fields. This is 
less than the nominal pixel scale (l'.'2) of the IRAC, ensuring 
the almost perfect pointing accuracy of the Spitzer observa- 
tions. Combining all uncertainties, a conservative estimate 
of the la referencing uncertainty is < Cf.'4 in both coordi- 
nates for all IRAC images in question. This is comparable 
with astrometric uncertainties of available Chandra/ ACIS 
and HRC images of the same fields in X-rays, but about 
twice as poor as for optical and near-IR images obtained 
with large ground based telescopes and the HST. Neverthe- 
less, this allows us to identify positionally on a subarcsec- 
ond accuracy level the mid-IR, near-IR, optical and X-rays 
objects in the respective pulsar frames, which have been 
obtained with significantly different spatial resolutions and 
pixel scales. 



3 RESULTS 

3.1 Morphology of the Vela field 

3.1.1 Identification of the Vela pulsar 

In Fig. [1] we compare fragments of the 3.6 and 5.8 /xm 
Spitzer images containing the Vela pulsar with the J 
and H band images of the same field obtained with the 
VLT/ ISAAC in the period of December 2000 - January 
2001 (|Shibanov et alj 120031 ). The point-like near-IR pulsar 
counterpart is marked by a thick line, and it is clearly vis- 
ible near the same position in the mid-IR. The measured 
coordinates of the counterpart candidate in the 3.6 fim im- 
age, where it is detected with a higher significance than at 
5.8 fim, are RA=08 h 35 m 20 S .635±0 S .042 and Dec=-45° 
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Figure 1. The ~ 20"x25" fragments of the infrared images of the Vela field obtained with the VLT/ISAAC in the J and H bands (left 
column), and with the Spitzer /IRAC (AOR 11542784) in 3.6 and 5.8 /an bands (right column). The Vela pulsar counterpart is marked 
by a thick line and the arrow shows the direction of its proper motion. Crosses are the expected pulsar radio positions defined separately 
for the VLT and Spitzer observation epochs. A knot-like object 'ol', an extended feature SE of the pulsar, possible associated with the 
X-ray PWN jet and counter-jet, and two background stars 'o2' and 'o3' are notified in the H band image. The contours in the 3.6 fim 
images are from the J band image. Henceforth the image pixel value scale is linear. 



10' 35'.'48±0'.'45 (J2000). The errors account for w 0'.'l2 ob- 
ject position uncertainties in the image and the mid-IR as- 
trometric referencing uncertainties estimated in Sect. 2.2. 
Within la the coordinates are consistent with the position 
of the near-IR counterpart, RA=08 h 35 m 20 S .628±0 S .01 and 
Dec=-45° 10' 35''16±07l (J2000), which is defined with a 
higher accuracy due to a better spatial resolution of the 
VLT observations obtained with the pixel scale of 0'.'147. 
The contours overlaid on the 3.6 /im image from the J band 
image underline the positional consistency of the detected 
mid-IR candidate and the near-IR counterpart. They show 
also that all background near-IR sources, except of the two 
faint near-IR point-like objects 'ol' and 'o2' located in im- 
mediate vicinity north-west of the pulsar, are identified in 
the mid-IR and vice versa. Hence, the appearance of a new 
mid-IR background source at the pulsar position, not vis- 
ible in the adjacent near-IR band, is very improbable. An 
apparent emission region extended from north-east to south- 
west in ~ 8 arcsec south-east of the pulsar is the stray light 
artefact from a bright saturated star located in ~ 0l5 north- 
east of the presented frames. This has been mentioned in 



Sect. 2.1 and may somewhat increase the backgrounds level 
near the pulsar position and decrease the significance of the 
candidate detection. 

The pulsar proper motion i s /x = 58±0.1 mas y r~ with 
a positional angle of 301?0±1?8 (|Dodson et all2003l ). Hence, 
the pulsar shift during four years between the VLT and 
Spitzer observation epochs is ~ / /23, that is indistinguish- 
able within the derived positional error budget. Attempting 
to improve the accuracy of the positional identification we 
have performed the relative astrometry. The VLT H band 
image with a point-source FWHM as C/.'5 was used as a refer- 
ence for finding the plate solution for the 3.6 fim image with 
FWHM w 1'.'8. The IRAF geomap/geoxytran tool with a 
general geometry option, allowing for the image shifting, ro- 
tation, pixel scaling, and a skew, was applied. A dozen of un- 
saturated isolated stars, which are common for both frames 
and are located in vicinity of the pulsar, were selected as 
reference points. The formal rms uncertainties of the plate 
solution fit were < (/.'09 with maximal residuals of < Cf.'l7 for 
both image coordinates. The pixel scale ratio was « 0.245, 
that coincides with the nominal one for the two images. Us- 
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ing this solution we transformed the 3.6 /im image coordi- 
nates of the mid-IR counterpart candidate to the H image. 
As a result, we obtained that the candidate counterpart is 
apparently displaced by 0'.'41±0 / . / 21 to the south from the 
near-IR pulsar counterpart. Similar displacement can be no- 
ticed also from absolute astrometry results, which have been 
used in Fig. [1] where the mid-IR candidate position in the 
3.6 /im image is also slightly shifted toward the south from 
the near-IR pulsar contour. In any case, the shift significance 
is too small, < 2a, to consider it seriously as a real one. 
Nevertheless, this may be a signature of a systematic shift 
caused by a faint extended feature located in the immediate 
vicinity southward of the pulsar and marked as a 'counter- 
jet' in the H band panel of Fig. [1] The feature is likely 
associated with the pulsar PWN and is clearly detected in 
the near- IR |Shibanov et alf2 003). but it cannot be resolved 
from the point source in the mid-IR due to a poor Spitzer 
spatial resolution. We shall discuss the implications of that 
later in Sect. 4.1.2. Using the plate solution and absolute 
astrometry of the H band image, we obtain the candidate 
coordinates, RA=08' 1 35 m 20 s .626±0 a .025 and Dec=-45° 
10' 35'.'51±(y.'26 (J2000), which are consistent with those 
obtained from the direct mid-IR absolute astrometry, while 
they are formally about twice more accurate than the latter 
ones. However, these improvements do not still allow us to 
distinguish the counterpart proper motion. The new formal 
position error budget is still comparable with the proper mo- 
tion shift of the pulsar, and including the possible systematic 
shift discussed above increases the positional uncertainties 
to the values obtained early from the absolute astrometry. 

The expected radio coordinates of the pulsar at the mid- 
IR observation epoch are RA=08 h 35 m 20 s .588±0 S .038 and 
Dec=-45° 10' 34'.'727±0'.'4 (J2000). They are derived from 
the most accur ate VLBI radio pos ition and proper motion 
measurements l|Dodson et al.ll2003h , and the errors are dom- 
inated by the mid-IR astrometric referencing uncertainties. 
The displacement between the expected position and the 
mid-IR candidate centre is ~ 0'.'9, or about (2-3)cr of the 
coordinate uncertainty level. This can be resolved in Fig. [T] 
where radio positions are marked by crosses. The direction 
of the offset approximately coincides with the pulsar proper 
motion positional angle. Considering available optical/near- 
IR images obtained early with the VLT and HST and ref- 
erenced with the USNO catalogue we got similar results: 
the expected radio positions are always marginally shifted 
by about 0'.'2-0'.'4 ahead of the counterpart positions ap- 
proximately along the line of the pulsar proper motion (cf. 
Fig. QJ. The variation of the offset values can be explained 
by different spatial resolution of the observations, while the 
same direction of the offsets is noticeable. At the same time, 
the relative values, such as the counterpart motion and its 
shifts from one epoch to others, are compatible with the ra- 
dio proper motion. Dodson et al. (2003) have discussed that 
some differences between the radio and optical absolute as- 
trometries of the pulsar can be due to different astrometric 
standards used in the two ranges, i.e., very accurately es- 
tablished extragalactic standards in the radio and mainly 
Galactic standards in the optical, where proper motions are 
sometimes poorly known. Therefore, the apparent offsets of 
the optical/IR counterpart from the radio position of the 
pulsar are likely caused by unknown systematic uncertainty 
of the optical astrometric standards for the Vela field. 




21.0 8:35:20.0 19.0 

RA I las 

Figure 2. The ~ 50"x85" fragments of the star-subtracted 
3.6 fim (top) and the X-ray (bottom) images obtained with the 
Spitzer /IRAC (AOR 11542784) and Chandra /HRC-l, containing 
the Vela pulsar and the brightest part of its torus-like PWN. 
The circle marks the pulsar position. The contour indicating the 
position of the two bright arcs and south-east jet of the X-ray 
PWN at the bottom panel is overlaid on the mid-IR image. 



Summarizing this part, the optical/near-IR counterpart 
of the Vela pulsar is firmly established based on its peculiar 
spectrum, proper motion, and pulsations with the pulsar 
period. Therefore, the la agreement of the mid-IR counter- 
part candidate position with the near-IR counterpart coor- 
dinates, obtained by us using both the relative and absolute 
astrometrief], allows us to ignore the radio-optical standard 
systematics discussed above, and to conclude that the source 
detected in the two IRAC bands is likely to be associated 
with the pulsar. To confirm the identification by the mid-IR 
source proper motion measurements, higher spatial resolu- 
tion and/or new Spitzer mid-IR observations at later epochs 
are necessary. 



3.1.2 Searching for the Vela PWN 
Immediate pulsar vicinity. 

Previous attempts to find any o ptical counterpart of the Vela 
PWN have been unsuccessful (|Mignani et aljfeoost ). In the 
near-IR H band image of Fig. [1] one can resolve a nearby 
knot-like source 'ol' north-west of the pulsar and a faint 
structure, 'counte r-jet', extended by a f ew arcsec just south- 
east of the pulsar (|Shibanov et alj|2003l ). Both are projected 

4 We got the same results considering another available optical 
images of different epochs as reference frames. 
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Figur e 3. The large-scale structure of the Vela plerion at 5 GHz, 8 fim, and in soft X-rays, as seen with the ATCA (right llDodson et al.l 
l2003al) '), Spitzer/IRAC (AOR 11374848) (Jeft), and Chandra /HRC-I (middle). The contours in the 8 /im image are from the 5 GHz and 
X-ray images. The central black and white contours in the X-ray and radio/mid-IR images indicate the location of the internal torus-like 
structure of the PWN zoomed in Fig. 2. Ticks in the 5 GHz image show magnetic field lines. 



on the pulsar X-ray jet/counter-jet structures, and may be 
associated with the PWN. We do not detect them in the 
mid-IR. A possible reason is that the second structure can 
be hidden in the pulsar-candidate source wings due to a 
lower spatial resolution of the Spitzer/IRAC as compared 
to the VLT/ISAAC (see Sect. 3.1.1). The more distant knot 
could be resolved, however it is not detected. It is also not 
visible in any available optical images obtained at different 
epochs. If it is not an artefact, it may be a time variable 
source. 

To search for PWN structures at larger scales, we 
performed PSF subtraction of background stars from the 
Spitzer images making use of the IRAF daophot/allstar 
tool. The results for the brightest inner torus-like part of 
the PWN are shown in Fig. [2] where we compare the star- 
subtracted 3.6 /im image with the archival Chandra /HRC-I 
X-ray imag^f]. The pulsar counterpart candidate was not 
subtracted and its position is compatible with the X-ray po- 
sition of the pulsar marked by the circle. Some remains of 
a poorly subtracted saturated field stars are visible, but no 
counterparts of the X-ray torus-like structure with jets and 
arcs are detected. We have not found them in other Spitzer 
bands as well. 

PJerion. 

The inner torus-like PWN, shown in Fig. [2] is only a tiny 
part of the Vela plerion, that extends in the radio and X- 
rays to several arcminutes. The respective field is shown in 
Fig. [3] where we compare the 8 /jm, soft X-ray, and 5 GHz 
images. We r e produc ed the 5 GHz image from fig. 10 of 
iDodson et alj (|2003al ). To reveal the plerion X-ray struc- 
ture the Chandra /HRC-I image is binned by 16 pixels and 
smoothed with 3 pixel Gaussian kernel. The image inten- 
sity scale is also different, as compared to that of Fig. [2 to 
reveal fainter emission from outer PWN regions, while its 
inner torus-like part is indicated by the contour taken from 
Fig. [5] to demonstrate its small extent as compared to the 
whole plerion size. We have applied contours from the 5 GHz 
and X-ray images on the 8 image. Contour taken from 



5 Obs. ID 1996, date 2002-01, exposure 49.5 ks, PI D. Helfand. 



the X-ray image shows the outer boundary of the X-ray ple- 
rion, where it merges with backgrounds, as it is seen with 
the HRC. Contour taken from the 5 GHz image indicate the 
outer boundary of the radio plerion, which consists of two 
bright, the north and south, radio 'lobes', located symmet- 
rically with respect to the inner PWN symmetry axis. The 
inner PWN X-ray contours are overlaid on the 5 GHz image, 
though this part of the PWN is not detected in the radio 
range. In contrast to the torus-like structure, whose sym- 
metry axis is directed north-west approximately along the 
proper motion of the pulsar, the plerion is mainly extended 
south-west. 

As well as in the radio image, the pulsar and its inner 
PWN are not detected at 8 ^m. However, the outermost east 
boundary of the large scale plerion and its south-west exten- 
sion in the radio and X-rays likely correlate with the west 
boundary of a bright extended emission structure in the mid- 
IR, which also extends south-west and, possible, originates 
from the dust formed from the remnant ejecta and heated 
by the PWN. In other SNRs at longer wavelengths the radi- 
ation from such extended structures is typically dominated 
by continuum emission from the dust, that can significantly 
contribut e to the radiatio n of a remnant at 8 /jm as well 
(see, e.g.. lRho etalll2009l ). Indeed, the similar but brighter 
structure is visible at 24 and 70 fim images obtained with 
the MIPS. However, they have much smaller and not rep- 
resentative FOVs, and we do not show them here. Never- 
theless, the fact, that the extended mid-IR emission is vis- 
ible in the images obtained with different Spitzer detectors 
in three different bands, at different AORs and observation 
conditions, suggests that it is real, but not a result of an 
artificial background variation. At shorter wavelengths the 
extended emission disappears, as expected for a rather cold 
dust. 

The centre of t he associated with the p lerion TeV source 
HESS J0855-455 l|Aharonian et all l2006h is located tens 
arcminutes south-west and is outside the frames shown in 
Fig. [3] It is mainly extended in the same south-west direc- 
tion as the X-ray/radio plerion. Unfortunately, there are no 
Spitzer images to combine a complete large scale picture 
of the extended mid-IR emission, around the Vela plerion 
including the TeV source. Nevertheless, a remarkable coin- 
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Figure 4. Top panels: An ~ l.'2xlf2 overview of the Geminga field, as seen with the Spitzer (AORs 19037696+19037952) at 3.6 /im 
(left) and Subaru in the / band (right). The / band image is smoothed with four pixel Gaussian kernel, that is about the mid-IR pixel 
size. It is partially affected by several vertical stray light artefacts form oversaturated stars. Circles mark the Geminga optical counterpart 
and its mid-IR counterpart candidate. The ~ 15" X 20" box region of the pulsar neighbourhood is enlarged in Bottom panels, where the 
Geminga counterpart is indicated by a line, and the arrow shows the direction of its proper motion. The contours in the bottom 3.6 ^m 
image, where the only AOR 19037952 is implemented, are from the / band image. 



cidence of its complicated east X-ray and radio boundary 
with the west boundary of the extended mid-IR emission 
suggests, that we likely see an interface between the pulsar 
wind and the SNR ejecta supported by a pressure balance 
between these two different fractions of the remnant. It is 
most clearly visible in the north-east part of the plerion, 
where a bright mid-IR protrusion fills a deep nebula bay 
between its north and south radio lobes. As seen from the 
radio image of Fig. [31 the magnetic field lines also lie along 
the protrusion/bay boundary, possible contributing to the 
pressure balance in this region. We see no extended mid- 
IR emission along the west plerion boundary. However, the 
magnetic field lines in the 5 GHz image are concentrated 
towards this boundary and mainly directed along it. The 
magnetic fields are stronger there and may play a key role 
in the shape and stability of the west part of the plerion. 



3.2 Morphology of the Geminga field 

3.2.1 Possible identification of the Geminga pulsar 

In Fig. [4] we compare fragments of the Geminga field im- 
ages obtained in the optical / band on January 2001 with 



the Subaru telescope and at 3.6 p with the Sp i tzer /I R AC. 
The / band image is taken from Shibanov et al.l ((2006) . The 
pulsar optical counterpart is indicated by the yellow circle in 
the top panels and by the thin line in the bottom ones, where 
an overview of the field and pulsar immediate vicinities are 
shown, respectively. The direction of the pulsar proper mo- 
tion of 178.2±1.8 ma s yr _1 with a position angle of 52?9±0?4 
|Fahertv et alj|2007h is shown by the arrow in the bottom 
I band image. For the 3.6 /im observation epoch the ex- 
pected displacement of the pulsar from its / band position 
is ~ l'.'22. It is consistent at la astrometric uncertainty level 
((f.'A) with the positions of a faint mid-IR source marginally 
detected at 3.6 fim, that is also marked as 'psr'. The position 
of the source itself is defined with a twice worse accuracy of 
~ 0"8, as compared to the Vela counterpart candidate case. 
The main reason is its faintness and a poor approximation 
of the object brightness distribution by a point-like source 
profile. Nevertheless, it is fully compatible with the pulsar 
proper motion path, and can be considered as a likely mid- 
IR counterpart of Geminga. 

Despite a worse spatial resolution of the Spitzer com- 
pared to the Subaru, as in case of the Vela pulsar, practically 
all background optical objects in the frames shown in Fig. [4] 
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Figure 5. The large scale overview of the Geminga field at 5.8 fim (left) and in X-rays (right). The contours in mid-infrared image are 
from the X-ray image. Bright components of the clumpy Geminga tail, as well as the two background IR/optical sources, 'si' and 's2', 
are marked. Geminga, marked by the cross in the X-ray image, is not detected in this infrared band. A bright X-ray blob B and other 
less bright blobs have no infrared counterparts. While the blob C at the end of the X-ray tail is likely associated with a bright star, and 
the blob Y may be associated with a red background source not visible at shorter IR/optical wavelengths. 



can be identified in the mid-IR and vice versa. This is un- 
derlined by black contours from the / band image, which are 
overlaid on the bottom mid-IR image. We have not found 
also any additional background object in immediate vicin- 
ity of Geminga in the near-IR image of th e field obtained 
with the HST in the adjacent F160W band (|Shibanov et all 
2006). Considering much larger FOV, as in the top panels 
of Fig. |4l we found only a few new faint mid-IR sources, 
which are not visible in the / band and located far away 
of the pulsar. Hence, the appearance of a new background 
source at the expected Geminga position in the mid-IR is 
very unlikely. This supports the correct identification of the 
pulsar. 

Unfortunately, the suggested counterpart candidate, 
is detected only at ~ 2<r level and only in the one 
Spitzer band. Combining the data from the two AORs 
(19037696+19037952), as is in the top panel of Fig. H has 
not allowed us to increase the significance of its detection 
due to a high Zodiacal background in the first AOR, though 
the source is seen at the same coordinates in both AORs 
(see Sect. 2.1). Therefore, this source has to be considered 
with a caution and needs a confirmation by deeper obser- 
vations. Nevertheless, summarizing all pros and contras, we 
will assume below, that the source is the real counterpart of 
Geminga, and estimate how this may affect the broadband 
spectrum of the pulsar. 

3.2.2 Searching for the Geminga PWN 

Deep X-ray observations have allowed to find fragments of a 
bow shock and an axial tail of Geming a at the distances of 
several arcminutes behind the pulsar (|Caraveo et al.l [20031 ; 
iPavlov et ai1l2010l ). We have not found any significant mid- 
IR counterpart of the bow shock at these scales. At the same 
time, the tail was found to consist of a time- variable clumps 
appea ring and disappear ing in various places from time to 
time (|Pavlov et al.ll2010h . and a few possible counterparts 
of the clumps can be considered. An overview of the pulsar 



field containing the whole extent of the axial tail is shown in 
Fig. [5] where we compare the mid-IR 5.8 /im and ACIS-flX- 
ray images. To go as deep as possible, the mid-IR image was 
combined from the two longest observations listed in Table 
1 (AORs 19037696 and 19037952). The brightest fragments 
of the X-ray tail detected a t the ACIS-I o bserv ation epoch 
are marked by B and C, in IPavlov et alj l|2010h notations, 
and yet another one by Y. The X-ray contours are overlaid 
on the mid-IR image. 

There is no any significant mid-IR counterpart for the 
X-ray time-variable B substructure. It was not seen in the 
early ACIS-S^ X-ray image and, hence, has been born at a 
time within a 3.5 yr interval between the two X-ray obser- 
vations. If we assume that appearing a bright X-ray clamp 
is accompanied with an optical/IR emission, then the ab- 
sence of B in the mid-IR suggests, that it was created either 
after the Spitzer observations and its age is < 9 months, or 
it has flashed before these observations, but the lifetime of 
the particles responsible for possible mid-IR afterglow asso- 
ciated with the X-ray flash is much shorter than that for the 
particles radiating in X-rays. 

The blob C nicely coincid es with a backgrou nd stellar 
source, as has been noticed by IPavlov et al l l|2010l ). We can 
add, that in the mid-IR this source is not a single but likely 
a bland of one bright an 2-3 fainter stars. 

The blob Y positionally coincides with a faint back- 
ground source in the 5.8 ^m image. The source is not visible 
at shorter wavelengths and hardly visible at 8 /im. This blob 
is close to a bright blob A, that was clearly visible in the 
early ACIS-S image (see, IPavlov et ai1l2010l ) and not seen 
in the ACIS-I image shown here, as well as in any available 
IR/optical image. It is not clear, whether this is a remnant 
of A, that has faded and moved by ~ 3" to the east from its 
initial X-ray position, or it is an independent substructure, 



6 Obs. ID 7592, date 2007-08, exposure 72.1 ks, PI G. Pavlov. 

7 Obs. ID 4674, date 2004-02, exposure 18.8 ks, PI D. Sanwal. 
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Table 2. The results of the aperture photometry for the pre- 
sumed mid-IR counterparts of the Vela and Geminga pulsars. 
The 3rd column shows the magnitudes and upper limits mea- 
sured using optimal aperture radii presented in the 2nd column. 
The presented aperture correction factors are for the flux. 
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that was flashed later in the time interval between the two 
X-ray observations. Additional X-ray and mid-IR observa- 
tions are necessary to understand whether the Y blob and 
the suggested IR counterpart are variable and have the same 
origin or not. If we indeed detected this blob almost simul- 
taneously in the X-ray and mid-IR, it could be younger than 
the blob B. The rest fainter and less extended X-ray sources 
possible related to the extended tail are not identified in the 
mid-IR. 

Consider i ng im mediate vicinity of the pulsar, 
iPavlov et ail (|2010l ) discussed a bow shock structure 
of the X-ray PWN with possible 'forward-jet', which is 
directed north-east of Geminga along its proper motion. 
As seen from Fig. O the forward-jet substructure, marked 
as 'si', is likely not related to the PWN, but is the X-ray 
counterpart of a red star-like background optical object, 
that is also detected in the near-IR and optical images 
(cf. Fig. [J) . In the deep ACIS-I image its profile overlaps 
with the outermost region of the nebula producing a false 
impression on the presence of a jet. The object si is also 
visible as a point-like isolated object in the ACIS-S image, 
where the exposure was short. A similar situation is for an 
apparent south-east 'protrusion' seen only in the ACIS-I 
image and marked as 's2'. It is likely to be a background 
source as well, that is brighter in the optical/IR, but fainter 
in X-rays as compared to 'si'. 



3.3 Photometry of counterpart candidates 



aperture photometry was performed on the star-subtracted 
IRAC images. 

For magnitude estimates we used IRAF daophot tasks 
with zero points given in the IRAC Instrument Handbook 
Prior photometric measurements, we converted each image 
from MJy steradian" 1 to DN units using the conversion fac- 
tors from image fits headers. In this case the daophot mag- 
nitude error output can be directly used for the flux uncer- 
tainty estimates. To find an optimal aperture radius based 
on the maximal signal to noise ratio, photometric growth 
curves were obtained and used. For the pulsar candidates at 
3.6 and 5.8 /im the optimal circular aperture radius was 4 
pixels (at 0'.'6 pixel scale) with the background annulus and 
dannulus in the ranges of 5-6 and 4-6 pixels, respectively. 
The measured fluxes were multiplied by aperture correction 
factors taken from the IRAC Data Handbook We have 
performed the PSF photometry as well and got the results 
consistent with the aperture photometry at la level. 

Conservative 3<r point source magnitude upper limits 
can be placed for the counterpart candidates in the mid- 
IR bands where the pulsar fields have been observed but 
no candidates were found. To do that, the circular aperture 
radii were selected encapsulating > 80 % of a potential point 
source flux, the apertures were centred at the pulsar posi- 
tions, and standard deviations and inherent flux variations 
within these apertures were used. The magnitudes were con- 
verted into fluxes in physical units and the results are sum- 
marized in Table(2] For Geminga at 3.6 /xm, where the coun- 
terpart candidate is visible only at 2a level, the data from 
the AORs 19037952 and 19037696 (Table [TJ provide consis- 
tent results. However, the flux uncertainty is smaller for the 
former AOR, where the Zodiacal background is lower, and 
we use this result for our further analysis. 



3.4 Multiwavelength spectra 

Let us assume that our mid-IR identifications are real, and 
consider, how this can affect the spectra of the pulsars. Using 
the measured mid-IR fluxes, available optical and UV data, 
in Fig. [6] we compile dereddened multiwavelength spectra 
of the Vela and Geminga pulsars in the mid-IR-UV range, 
and compare them with the respective spec t rum o f the Crab 
pulsar taken from ISandberg fc Sollermanl l|2009l ). The in- 
terstellar extinctions E(B-V) are 0.0 55 and 0.023 mag to- 
ward s Vela and Geminga, respectively (|Shibanov et al ]|2003l . 
2006). These are small enough, and dereddening corrections 
in the mid-IR bands are negligible as compared to the mea- 
sured flux uncertainties. For instance, the magnitude cor- 
rection Aa at A = 3.6 fim is < 0.01 mag for both objects. 
For clearness, in Fig. [5] we notify the broad bands, where 
pulsar fluxes were measured, and indicate which ones were 
obtained with the Spitzer /IRAC or with the MIPS. 

In contrast to the Crab, that shows almost a flat spec- 
trum in the considered range, the compiled spectrum of 
the Vela pulsar demonstrates a strong, an order of mag- 
nitude, mid-IR flux excess, that progressively increases with 
the wavelength, as compared to the optical range. Geminga 



To minimize the flux contamination of the detected pulsar 
counterpart candidates from nearby background objects the 



8 see, http://ssc.spitzer.caltech.edu/irac 

9 IRAC Data Handbook, table 5.7. 
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tive ones. We consider this as an additional evidence, that 
the detected mid-IR emission is indeed associated with the 
pulsars. 

Accepting that, we can now update the rotational phase 
averaged spectral pictures for both pulsars in the whole 
observational range, from the radio to hard 7-rays (cf. 
Shib anov et alj|200al2006h , including our mid-IR data, and 
the data obtained r ecently with the Fermi 7-ray observatory 
jAbdo et al.ll2010al lbh. 

As seen from Fig. the flux density of the mid-IR emis- 
sion from the Vela pulsar becomes comparable with the den- 
sity of the spectral bump produced by thermal emission from 
the atmosphere (atm.) at the surface of the NS visible in soft 
X-rays. At the same time, within large uncertainties it may 
be compatible with the long wavelength extrapolation of the 
nonthermal soft X-ray spectral tail detected with Chandra 
and the RXTE and described by a power-law (hatched re- 
gion) . The mid-IR flux upper limits are still lower than flux 
densities of a coherent radiation at radio frequencies. 

Similar spectral behaviour is observed for Geminga. The 
main difference is that Geminga is fainter in the optical/IR, 
and the observed mid-IR excess is much less significant. Nev- 
ertheless, the steep IR flux increase towards low frequencies 
can not be excluded by the derived mid-IR flux upper lim- 
its. This suggests, that at several 10s microns the mid-IR 
flux density might be comparable with the soft X-ray ther- 
mal blackbody-like emission bump from the Geminga sur- 
face (BB) and compatible with the low frequency extension 
of the nonthermal soft-X-ray spectral tail, as in the Vela 
case. 

However, from the spectral behaviour of both pulsars 
shown in Fig. [6] and Fig. [JJ it is most likely that the Spitzer 
mid-IR observations reveal us a new spectral component in 
the emission of these rotation powered-pulsars, whose origin 
is different from those of the rest components detected early, 
and, perhaps, from that of the Crab in the mid-IR. 



Figure 6. Unabsorbed spectra of the Crab, Vela, and Geminga 
pulsars f rom the mid-IR to th e opti cal— UV. The Crab data 

~" teoosft. 



and near-IR data are from 


Romani et alj J2005I1. iMignani et al.l 


1120071), and IShibanov et al 


d2003t), respectively. The Geminga 



UV and optical/near-IR d ata are from Karealtsev et all (2005) 
and IShibanov et al.l (I2006T). respectively, the optical spectrum is 
from lMartin. Halpern &: Schiminovicbl lll998h . 



shows a similar, though much less significant, excess, as fol- 
lows from its marginal detection. It is important to note, 
that an intentio n of the flux inc r ease t oward s the IR has 
been noticed by IShibanov et al. 

I (|2003l . l200fj ), when both 
pulsars were identified in the near-IR JH bands. However, 
the significance of that was low, and the presence of such a 
strong mid-IR excess was difficult to expect, based on the 
flat spectrum of the Crab pulsar, that has been so far the 
only rotation-powered pulsar detected in the mid-IR. The 
mid-IR fluxes of the Vela and Geminga candidates are likely 
to be in agreement with their near-IR counterpart fluxes, 
confirming the reality of the apparent flux increases towards 
the IR, resulted early only from the near-IR data. Thus, 
the compiled spectra appear to be real, but not the tenta- 



4 DISCUSSION 

Subarcsecond positional coincidence of the detected point- 
like object at 3.6 and 5.8 /im with the Vela pulsar coor- 
dinates suggests, that this object is likely to be the mid-IR 
counterpart of the pulsar. This is also supported by the con- 
sideration of its spectral energy distribution (SED), showing 
that the measured mid-IR fluxes are compatible with the 
long wavelength spectral excess observed early in the near- 
IR. The same is likely to be true for fainter Geminga, with 
a caution that its mid-IR counterpart candidate is detected 
only marginally and solely in the 3.6 /im band. Additional 
confirmations can be obtained from near-IR detections of 
both pulsars in the K band, where they have not been ob- 
served, and by deeper mid-IR observations, to detect the 
Geminga counterpart with a higher significance and to es- 
tablish both identifications by proper motion measurements 
of the counterparts in this range and checking their consis- 
tence with the proper motions of the pulsars. Extension of 
the pulsar observations to the far-IR would be also useful in 
this context. 
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Figure 7. Multiwavelength unabsorbed spectra for the Vela and Geminga pulsars from the radio to hard 7-rays compiled from the data 
obtained with different instruments, as indicated in the plots. Mid-fR flu xes are marked by red circ les. Down triangles mark upper limits. 
Filled squares in the radio range of Geminga are its radio detections bv lMalofeev fe Malovl il997T l. 
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4.1 The nature of the mid-IR excess 

If our mid-IR identifications of the two pulsars are real, the 
question is what is the nature of the observed mid-IR flux 
excesses. 

Besid es the Crab, where the ex cess is likely to be 
marginal dSandberg fc Sollermanl [20091 there are only two 
other isolated NSs d etected in the mi d-IR. These are the 
AXPs 4U 0142+61 (|Wang et all 120061 ) and IE 2259+586 
( Kapl an et al J 12009! ) . Both have been widely considered as 
radio-silenll 10 !. and belong to a class of magnetars, i.e., NSs 
with extremely strong magnetic fields. They have significant 
mid-IR excesses, which have been interpreted as the emis- 
sion from X-ray irradiated fall-back discs around the NSs. 
Such discs can be formed aro und isolated NS s through a 
process of supernova fall-back (|Chevalierjll989l ) . 

The marginal excess over the power-law extrapolation 
of the optical SED for the Crab was suggested to be caused 
by a nearby PWN knot located in 0'.'6 from the pulsar. The 
knot is much fainter than the pulsar in the optical, but it 
has a very red spectrum with a steep flux increase towards 
the mid-IR, where the knot brightness can be comp arable 
to that of the pulsar l|Sandberg fc Sollermanl [2009). This 
can result in an overestimation of the measured pulsar flux 
in this range, since the Spitzer spatial resolution does not 
allow one to resolve the knot from the pulsar. 

A planet, or a cool low-mass stellar companion could 
also provide an excess. However, this is unlikely for our pul- 
sars, since no planets or binary companions have been re- 
ported from timing observations of Vela and GemingJH. 

Finally, the observed excesses can be produced in mag- 
netospheres of the pulsars. The latter possibilit y has also not 
been excluded for the magnetar IE 2259+586 ( Kaplan et al.l 
2009). Below we discuss in more detail the possibilities listed 
above for the interpretation of the Vela and Geminga SEDs. 
We will focus mainly on the Vela pulsar, whose mid-IR coun- 
terpart candidate is firmly detected. 



4-1.1 The fall-back disc 

As protopl anetary discs around ordinary pre-main sequence 
stars (e.g., iDullemond fc Monnierl |2010T) or dust discs de- 
tecte d around white dwarfs (e.g., iMelis et"alll201ol ; iFarihil 
2010), the fall-back discs around pulsars can lead to emission 
excesse s in pulsar spectra in the mid-IR and submillimetre 
band s ([Foster fc Fisc her 1996; |p"erna. Hernquist fc Naravanl 
l2000l) . 

To find out whether the Vela and Geminga spectral ex- 
cesses can be explained by radi ation from hypothe tical discs 
around these objects, following IWang et al.l (|2006f ). we con- 
sider passive dust discs heated by X-ray irradiation from 
the pulsar s. We adopt a mode l of the irradiated disc emis- 
sion from IVrtilek et all ([l990h . It is assumed that there is 
enough gas to provide the pressure support for the dust, 
but its contributions to the heating balance due to a viscous 
dissipation and the emission are negligible. In this model the 
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Figure 8. The spectral fit of the dereddened Vela pulsar optical— 
IR SED by a sum of the fixed power-law component (dashed 
line), which describe contribution from the pulsar itself, and X- 
ray irradiated passive dust disc model (dotted line). The solid 
line is the best-fitting result and the hatched region shows ltr 
uncertainty of the fit. The dotted line is the contribution from 
the thermally emitting disc with the best-fitting parameters. For 
the parameters see the text. 



disc is optically thick and at a fixed disc radius it radiates 
as blackb ody with a temper ature depending on the radius 
oc R~ 3 / 7 (jVrtilek et allll990T ). The parameters determining 
the emission properties of the disc are the inclination angle 
f , albedo r\, and the internal and external disc radii Ri„ and 
Rex- All of them cannot be certainly constrained from the 
spectral fit to the observed IR SEDs containing only four 
points and upper limits. One has to select most critical ones 
to be found from the fit, and fix the others at reasonable 
values. Albedo strongly affect the results and in our fits we 
fixed it using a set of values between 0.5 and 1. We fixed 
also possible contributions from the pulsars in the IR, as- 
suming that in this range their spectra are consistent with 
a power-law extrapolation from the optical-UV range. The 
observable disc properties depend also on the X-ray lumi- 
nosity of the pulsar and its distance, which are known for 
our puls ars from previous s tudies: L\ ela = (5.3+0.5) x 10 32 
erg s" 1 dPavlov et all 1200 ll ). L% em = (3+1) xlO 31 erg s" 1 
jveza _ 00-7+19 pc (|bodson et all 



l|Kargaltsev et al.ll2005h . ~d^ ela = 287; 
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by recent analysis. 



l2003h . and d Gem = 250±£jj° pc (|Fahertv et al.ll2007h . 

For Vela the inclination angle can be fixed at 64°, ac- 
counting for the viewing angle of the pulsar rotation axis , 
estimated based on the PWN geometry (|Helfand et al 1 l200ll ; 
iNg fc Rom ani 2008). Within uncertainties this angle is com- 
patible with that der ived from the polarim etric radio obser- 
vations of the pulsar (| Johnston et al|[2005l ). We assume, that 
the disc and pulsar rotation axes coincide with each other. 
As a result, the irradiated disc model provides a good fit for 
the Vela SED at 77 < 0.85. An example is presented in Fig. [8] 
where we fixed 77 at 0.75. In this case Ri n = 0.28 ± 0.05 
R©, Rex = 1.4 ± 0.4 Rq, and the dust temperatures at the 
inner and external disc boundaries Ti n — 834 ± 63 K, and 
Tex = 418 ± 56 K, respectively. 

The Ri„ value derived from the fit is much larger than 
the pulsar light cylinder radius Rlc ~ 0.6xl0 -2 R0. This 
allows Vela to work as the radio pulsar. At the same time 
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it is about ten tim es smaller than that of AXP 4U 0142+61 
|Wang et alj|2006h . suggesting a more compact disc around 
the Vela pulsar than the disc around the AXP, while the 
inner disc temperatures for both NSs are similar and con- 
sistent with the sublimation temperature of dust > 1000 
K. This is not a surprise, since the X-ray luminosity of the 
AXP, 9xl0 35 erg s _1 , is about three orders of magnitude 
higher than that of the Vela pulsar, and the AXP dust disc 
can survive only at larger distances from the NS. Therefore, 
Rin derived from the fit is in agreement with the inner disc 
boundary, that can be estimated independently from the 
dust destruction te mperature, as it has also been noticed by 
IWang et all (|2006h for the 4U 0142+61 case. We also note, 
that our R ex value is compatible with a condition R ex < R© 
found to be c orrect for alm ost all dust discs observed around 
white dwarfs (|Farihill2010l ). This means that R ex cannot ex- 
ceed the Roche radius , within which any rocky object would 
be tidally disrupted (jMelis et al.ll2010n . The Roche radius 
for NSs, as for white dwarfs, is about Rg0 

If Rin 3> Rlc, then fall-back discs around isolated pul- 
sars can be swept away by radiation pres sure of the intensiv e 
magneto-dipole radiation of the pulsar I Shvartsmarj|l970j ). 
However, it was shown l|Eksi fc Alparl l2005l ~ that such a 
disc is not disrupted by this process, if Ri n < 2Rlc and 
Rin < 100J?ic in cases of the orthogonal and aligned mag- 
netic rotators, respectively. For Vela our best fit value Ri n ~ 
45Rlc lies in the middle of this range, providing stability of 
the disc against the radiation pressure at the angles between 
the magnetic and rotational axes of the pulsar < 10° . This 
is in contrast to a value of ~ 70° estimated f rom a combina- 
tion of the pulsar polarimetric observations Jjohnston et al. 
2005), and X-ray PWN geometry analysis dNg fc Romani 
20081 '). The latter value suggests the disruption of the disc, 



since at this angle we cannot obtain any reasonable fit to 
the Vela SED with Rin, that is close enough to the Rlc 
value to satisfy the stability criterion. Nevertheless, if the 
Vela pulsar is indeed closer to the orthogonal rotator, but 
not to the aligned one, the disc might still survive, if we 
assume that optically thin gas disc component exists inside 
Rin, as it is for dust discs ar ound pre-main-sequence stars 
( Dul lemond fc Monnierllioibl ). This component cannot con- 
tribute significantly to the IR emission, but may support the 
stability of the disc against the radiation pressure. 

At the same time, the problem of the dust disc structure 
and stability in the presence of a strong relativistic particle 
wind from energetic pulsars is still not resolved (see, e.g., 
iBrvden et al.ll2006l ; [jonesll2007l . and references therein). The 
presence of the p lanetary system orbiting PSR B1257+12 
(jWolszczanll 19941 ) strongly suggests that its formation must 
involve a preplanetary disc circling the NS. Despite of many 
efforts ijBrvden et alJl2006h . the most stringed upper lim- 
its on the dust disc emission in this system obtained with 
the Spitzer /MIPS, suggests that the dust, if it really exists 
there, is significantl y cooler than tha t derived by us for the 
Vela pulsar and bv lWang et all (|2006r i for AXP 4U 0142+61. 
Therefore, we conclude, that a simplified model of a passive 
dust disc can formally explain the IR excess for Vela, but 



The exact value depend on the densit y of a rocky object as 
well as on mass and radius of a star (Mc lis et al.ll20l"oh . 
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Figure 9. The comparison of the spectral energy distributions for 
the likely Vela pulsar counterpart and the two structures located 
in ~ 2" vicinity of the pulsar and projected on the NW X-ray 
PWN jet ('ol') and the SE counter-jet of the PWN (SE-structure 
), shown in the attached fragment of the H band image. The 
data-points of the structures are outlined by ellipses and notified 
in the top panel plot, while the rest points belong to the pulsar 
counterpart. Dashed lines show tentative spectral fits for these 
structures by power laws with similar spectral indices of about 
3.1. As seen, the spectral slope of the pulsar counterpart candidate 
in the IR is compatible with the spectral slopes of the structures. 
Large and small circles in the image indicate the FWHM of the 
PSF and the pulsar position with its Icr uncertainty in the mid- 
IR, respectively 



it is not clear, whether such a disc can survive around this 
energetic pulsar. 

We note also, that the Vela IR SED can be roughly 
fitted by a pure blackbody model with T ~ 597 K, and a 
blackbody radius of ~ 0.7 Kq. The latter one is too large 
for a planet or a cool brown dwarf. This is an additional 
argument against a cool binary companion of the pulsar as 
a source of the excess. 

We have carried out the similar analysis for Geminga. 
The results are very uncertain due to larger uncertainties of 
its IR fluxes. We can note only that the fall-back disc around 
Geminga cannot be excluded by current data. The inner disc 
radius seems to be comparable with the Geminga light cylin- 
der radius, Rlc ~ 0.016 Rq, perhaps, suggesting the disc 
stability, that m ay lead to unstable behav iour of Geminga as 
the radio pulsar ( Malofeev fc Malov 1997j). However, deeper 
IR studies are necessary to make any definite conclusions. 
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4-1.2 An unresolved PWN structure 

As in the Crab case, the IR excesses of Vela and Geminga 
may be produced by unresolved nearby structures of their 
PWNe. 

There are two faint red structures in ~ 2" vicinity of 
the Vela pulsar, whic h have been detected in the near-IR 
|Shibanov et al.l l2003h . They are projected on the X-ray 
PWN jet and counter-jet, and can be barely resolved in the 
bottom-left panel of Fig. [1] For clearness the respective re- 
gion is enlarged in the bottom panel of Fig. [9] A knot-like 
object 'ol' is located wl"2 NW of the pulsar and detected 
in the H band. An extended «l'.'5 structure is located just 
SE of the pulsar, partially overlapping with it, and detected 
in the J and H bands. Both objects are not seen in the opti- 
cal bands, and the SE structure has showed some signs of a 
short-time variability in the near- IR (|Shibanov et al.| [2003). 
The SED of 'ol' with Fj < 0.4 /ijy and F H = 0.97(18) pJy 
can be extrapolated to the mid-IR using a power-law with 
a spectral indej^l a u 3 (Fig. |9}. This fact, and the spa- 
tial separation of 'ol' from the puls ar, that is comparable 
to that of a nearby west star ('o3', in lShibanov et all (|2003l ) 
notations) visible in the mid-IR and optical (cf. Fig. [TJ, sug- 
gest that it could be detected in the Spitzer /IRAC bands 
at a similar brightness level as the pulsar candidate in the 
wing of its PSF, whose mid-IR position and FWHM are 
marked by circles in the bottom panel of Fig. [5] However, 
we do not detect it there. This means, that either 'ol' has 
a very peculiar and unrealistic SED with a sharp flux peak 
in the H band, or it is a time variable object. If 'ol' and 
the SE structure are associated with the Vela PWN, the 
latter assumption is compatible with a strong time variabi l- 
ity of the PWN observed in X-rays i|Pavlov et al.l l2001al ). 
As seen from Fig. [§] at the epoch of the near-IR observa- 
tions the integral brightness of the SE structure, with Fj = 
2.59(12) /ijy and F H = 4.94(38) /iJy, was higher than the 
brightness of the pulsar, while 'ol' was fainter. The spectral 
slopes of both structures are comparable to each other and 
to the spectral slope of the pulsar candidate in the mid-IR. 
Since the SE structure cannot be resolved from the pulsar 
in the Spitzer images, it can considerably contaminate the 
measured pulsar fluxes, assuming that the emission inten- 
sity of this structure in the mid-IR is compatible with the 
long wavelength extrapolation of its near-IR SED. This may 
be also supported by our astrometry, that shows a marginal 
shift of the mid-IR candidate position (black circle in the im- 
age of Fig. [9| from the near-IR counterpart centroid toward 
the SE structure. 

We can then speculate, that the SE structure can, at 
least partially, explain the observed excess, if this structure 
became fainter at the Spitzer observation epoch, as it is for 
'ol', but saved its steep spectral slope. In this case, we can 
expect a variation of the excess intensity with time, if the 
SE structure is indeed associated with the PWN, that varies 
with time. New Spitzer /IRAC observations in 3.6 and 5.8 
/xm, and high spatial resolution observations in the near-IR 
K s band would be useful to verify that. We note, that the 
spectra of 'ol' and the SE structure are apparently about 2- 
3 time s steeper than the spectrum of the Crab knot with ct v 
« 1.3 (|Sandberg fc Sollermanll2009h and the X-ray spectrum 



The index a„ is defined as F v oc v 



of the PWN. The re asons for that are not clear (see, e.g., 
IShibanovet~al]|2003l ). and further IR observations can help 
to understand their real nature. 

The object 'o2' is also not visible in the mid-IR, while 
it is detected in the optical RI and near-IR JH bands, 
showing that it is likely to be a distant (< 10 k pc) main 
sequence background star (|Shibanov et al.l l2003h . For in- 
stance, for MOV star with the observed colour J—H ~ 1 at 
a Galactic Ay ~ 4, we obtain d ~ 10 kpc and the expected 
fluxes at 3.6 and 5.8 fim of < 1.5-2.5 /xJy. This is bellow 
the Spitzer observation detection limits and about 20 times 
fainter than the candidate flux at 5.8 ^m. Therefore, the 
much brighter candidate can not be significantly contami- 
nated by 'o2', even if its PSF overlaps with this star. 



4-1.3 Pulsar magnetospheric emission 

As seen from Fig. [7] the nonthermal magnetospheric emis- 
sion dominates multiwavelength spectra of both pulsars 
from the radio to 7-rays, except of a narrow range in soft 
X-rays, where the thermal radiation from the surfaces of the 
NSs is more intensive. The spectra are mainly of a power- 
law, but their slopes are different in different spectral do- 
mains, and there are several spectral breaks, e.g., between 
the 7-rays and the hard X-rays, the soft X-rays and the 
optical, etc. These features are not yet clearly understood, 
and in this context, a new break between the optical and 
IR does not look very exclusive, and may represent an ad- 
ditional component of the NS magnetospheric emission. 

Except of the Crab, three, of four young PWNe so far 
identified in the optical-IR (PSR B0540-69, 3C 58, and 
G292.0+1.8. dWilliams et al.ll2008l;ISerafimovich et aLpOol 
Shibanov et al.ll2008l ; ISlane et al.ll200Sl ; IZharikov et al-lbood ; 
Zvuzin et al. 120091 ') show strong mid-IR excesses with steep 
power-law spectral slopes in the IR. It might be possible, 
that the relativistic particles responsible for these IR ex- 
cesses are ejected from magnetospheres of the pulsars pow- 
ering the PWNe. In this case, the same particles can provide 
the IR excesses in the magnetospheric emission of the pul- 
sars. 

However, only the detection of pulsations with pulsar 
periods would be a strong evidence on the NS magneto- 
spheric origin of an enhanced emission observed for Vela 
and tentatively for Geminga in the IR. Until the time, when 
timing observations of pulsars in the IR will be possible, 
imaging observation of other suitable rotation powered pul- 
sars in the IR would be useful, to expand the sample of NSs 
detected in this range and to verify, if the IR excess, like 
that of Vela and Geminga, is a more common feature in the 
pulsar emission, than a flat optical-IR SED of the Crab. 



4.2 Mid-IR PWN signatures 

The inner torus-like PWNe structures around the 
young Crab, J1124-5916, J0205+6449, and B0540-69 
pulsars are well identified in the optical and mid- 



IR (ISandberg & Sollermanl I200E 


; IZharikov et al.l 12008: 


Zvuzin et al. 20091: IShibanov et alj 


2008; Serafiniovich et al.l 


20041; 1 Williams et all l2008f). The latter three obiects are 



even brighter in the IR than in the optical, and their 
SEDs demonstrate a steep flux increase toward longer wave- 
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lengths. In all cases the optical/IR PWN counterparts have 
approximately the same sizes as in X-rays. 

The Vela and Geminga pulsars are considerably older 
and less energetic, as compared to the above group, and 
their X-ray PWNe are fainter as well. This is likely to be 
the main reason why we do not observe PWNe counterparts 
of these pulsars in the optical and IR at the same spatial 
scales as they are seen in X-rays. However, we cannot ex- 
clude the presence of a hint of a counter-jet fragment at a 
small distance from the Vela pulsar. As have been discussed 
above, such poorly resolved small scale structures can be a 
reason of the observed IR excesses in the pulsar SEDs. The 
compactness of these structures, if they are real, show that 
the inner PWN evolves with the pulsar age and fade more 
rapidly at longer wavelengths than in X-rays. We note also, 
that there are still no reports on the detection of the inner 
tours-like structure with jets for the Vela PWN in the radio 
range, while a 10 arcsec long Geminga tail s tructure was re- 
centl y detected with the VLA at 4.5 GHz (jPellizzoni et al.l 
|201lJ). The Geminga radio flux density is compatible with 
the long wavelengths extension of the tentative mid-IR ex- 
cess reported here. 

At the same time, the mid-IR observations allowed us 
to define more exactly the Geminga X-ray PWN shape in 
the immediate vicinity of the pulsar. Its forward-jet and the 
south-east protrusion are likely not to be the real PWN 
structures, but are X-ray counterparts of background red 
stellar object 'si' and 's2' (Fig. 0)- The same is true for the 
distant blob C, that coincides with a bright IR background 
object. More careful X-ray spectral analysis is necessary to 
confirm or reject the nature of the 'si' and 's2' objects, that 
is outside the scope of this paper. 

The comparison of the X-ray and 8 /im images of Vela 
shows a good coincidence of the west boundary of the ex- 
tended IR emission of the remnant with the east boundary 
of the large scale X-ray plerion (Fig. [3]). Further comparison 
of the mid-IR image with the high spatial resolution ra- 
dio images of the Vel a plerion, obtained with the ATCA by 
iDodson et al.l (|2003aF ). likely supports this statement. A bay 
between the north and south plerion lobes in the radio image 
is filled by a bright protrusion of the extended mid-IR emis- 
sion. The magnetic field lines lie along the bay /protrusion 
boundary, also suggesting that we likely see an interface be- 
tween the relativistic pulsar wind and a dust SNR ejecta 
provided by the pressure balance between these two differ- 
ent parts of th e remnant. No lobes structures are visible in 
the optical (cf., Mi gnani et al.|[2~003h and more careful anal- 
ysis is necessary to check this hypothesis. 

To summarize, we conclude, that Spitzer mid-IR obser- 
vations have opened a new spectral window for the study 
of the pulsars. The detections of the likely counterparts of 
the Vela and Geminga pulsars, and PWNe for several Crab- 
like pulsars reveal strong mid-IR excesses in the emission 
of these objects, that was rather difficult to expect in ad- 
vance. Further mid-IR studies of these and other pulsars are 
necessary to increase the number of the objects detected in 
the IR, and understand the nature of their emission in this 
range. 
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